Absence of non-linear Meissner effect in YBa 2 Cu 3 06.95 
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We present measurements the field and temperature dependence of the penetration depth (A) 
in high purity, untwinned single crystals of YBa2Cu3 06.95 in all three crystallographic directions. 
The temperature dependence of A is linear down to low temperatures, showing that our crystals are 
extremely clean. Both the magnitude and temperature dependence of the field dependent correction 
to A however, are considerably different from that predicted from the theory of the non-linear 
Meissner effect for a d-wave superconductor (Yip-Sauls theory). Our results suggest that the Yip- 
Sauls effect is either absent or is unobservably small in the Meissner state of YBa2Cu3 06.95- 
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PACS numbers: 74.72.Bk, 74.25. Nf 

The identification of a strong linear term in the tem- 
perature dependence of the magnetic penetration depth 
(A) of YBa2Cu3 06.95 was a key development in identify- 
ing the d-wave symmetry of the order parameter in high 
temperature superconductors. El The linear increase of A 
with increasing temperature in a d-wave superconductor 
reflects a decrease in supcrfluid density caused by thermal 
excitations near the nodes in the order parameter, and 
is insensitive to the direction of the superflow. Yip and 
SaulsaB have proposed that a much more sensitive test 
of the order parameter symmetry can be made by mea- 
suring the field dependence of A in the Meissner state. 
They argued that at sufficiently low temperature, a pure 
d-wave superconductor should have a weak, linear field 
dependence, the magnitude of which depends on the rel- 
ative orientations of the applied field and the nodes in 
the order parameter. 

Although the linear temperature dependence of A is 
now well established in many materials the field depen- 
dence has proved much more difficult to isolate. There 
is to date no conclusive evidence as to the presence or 
absence of this effect in the Meissner state in any ma- 
terial. In this paper we present measurements of both 
the field and temperature dependence of A below H c \ in 
high purity, untwinned single crystals of Yh^CusOg.gs, 
in all three crystallographic directions (H\\a,b,c). We 
find that in some crystals there is a relatively strong field 
dependence to A which might naively be attributed to the 
Yip-Sauls effect. However, a closer examination reveals 
that it is almost certainly of extrinsic origin. Our cleanest 
crystal shows a field dependence of A which is up to a fac- 
tor 7 lower than that expected from the Yip-Sauls theory 
and is not suppressed by increasing temperature. This 
suggests that the Yip-Sauls effect is either absent or un- 
observably small in the Meissner state of Yl^CusOg.gs. 

In the Yip-Sauls theory the field dependence of A at 
low temperature is a direct consequence of the presence 
of the nodes in the order parameter. In the presence of 
a finite superflow (v s ) the quasiparticle energy levels are 
Dopplcr shifted by an amount 8e oc v s -Vf, where Vf is the 
Fermi velocity. In a d-wave superconductor, levels near 



a node with vj in the opposite direction to v s will be 
shifted below the Fermi level and there will be a transfer 
of quasiparticles from comoving to countermoving states, 
thus resulting in a reduction of the effective superfluid 
density. In an s-wave superconductor the effect is absent 
at low temperature because the finite gap prevents the 
shifted levels from being occupied. Yip-Sauls calculate 
the field dependence in a d-wave superconductor to be 



AX(H) _ H 



(1) 



where, Hq is of order the thermodynamic critical field (Xu 
et alB estimated it to be 2.5T for YBa2Cu3 06.95), and a 
depends on the relative orientation of the nodes and the 
field. For field directed along an antinode, we calculate 
from this formula that the-slope dX a (H)/dH = 4.7 x 10~ 2 
A/Oe for YBa 2 Cu 3 6 .95,H or 4.7 A in a field of 100 Oe. 
The Yip-Sauls effect is very sensitive to temperature and 
sample purity. At finite temperature AA is only linear 
in H above a crossover field H* ~ ki,THo/2Ao which 
exceeds H c \ for Tj-fSj 3 K (here Ao is the maximum en- 
ergy gap - 1.9T C M Below H* AA - H 2 and dX/dH is 
much reduced. This high sensitivity on temperature can 
be used to distinguish the Yip-Sauls effect from other 
possible (non-intrinsic) contributions. 

Our measurements of A were performed by measuring 
sample induced changes in the inductance of a copper coil 
(containing the sample). The coil formed part of a tunnel 
diode driven LC oscillator operating at 13 MHz. Changes 
in the resonant frequency of this circuit are directly pro- 
portional (to within a calibration constant) to changes in 
the effective volume of the inductor. By careful design 
of this measurement .circuit (following the prescription 
given by Van Degrifta), we have achieved a RMS noise 
level in the oscillator of 1 part in 10 9 / y/Hz. For a typical 
crystal with 1mm 2 surface area, this translates to a noise 
level in the measurement of A of ~ 0.1 A. This resolutio, 
is about 100 times better than previous measurements 

A small superconducting solenoid coaxial with the 
probe coil provides a variable, static, dc field, up to 
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^1000 Oe collinear with the RF probe field. The com- 
plete experiment was placed inside a mumetal shielded 
Dewar. The remnant field inside this was measured, us- 
ing a fluxgate magnetometer, to be < 2 mOe. The RF 
probe field to which the sample is subjected is estimated 
to be also < 2 mOe. For a single crystal, with the H 
parallel to the plane, the total flux in the sample is of 
order one flux quantum ($o)- These very low field values 
help keep to an absolute minimum the number of vortices 
trapped in the sample as it is cooled through T c . 

The sample is located on a sapphire rod, the temper- 
ature of which is varied independent of the rest of the 
apparatus. The sample rod may be drawn out of the mea- 
surement coil at low temperature so that the background 
field dependence of the apparatus can be measured in- 
situ. We found that this background was highly repro- 
ducible at low temperature, but could change markedly 
upon thermal cycling, making an in-situ measurement vi- 
tal. There was a further background originating from the 
sapphire sample holder, which was temperature depen- 
dent, but highly reproducible between runs. This was 
therefore measured in a separate run with the sample 
removed. These backgrounds typically give a frequency 
shift corresponds to around a 3 A equivalent change in A 
for a mm sized single crystal in 100 Oe. 

The effective volume of the sense coil was calculated 
by measuring the shift in frequency when a small sphere 
of Aluminium was inserted. In the B\\a, b geometry 
the measured frequency shifts are then easily related to 
changes in penetration depth simply by measuring the 
area of the faces of the sample. This neglects the contri- 
bution of the currents flowing along the c-axis, although 
for the aspect ratio of the crystals measured here, this 
results in less than a 1% error in AA(T = 30K)R For the 
B\\c geometry a calibration factor is rather more difficult 
to calculate and so it was estimated by comparing the 
temperature dependence of the frequency shifts in this 
orientation to that for B\\a, b. This factor was then used 
to relate the field dependent shifts to AA(iJ). 

The single crystals measured were grown in Yttria sta- 
bilized Zirconia crucibles as described in Ref. ||. They 
were annealed so as to give optimal doping. The crys- 
tal used for the bulk of the measurements described here 
had dimensions a x b x c = 0.80 x 0.51 x 0.011 mm 3 , was 
naturally untwinned and had a T c of 91. 4K. The crystal 
contained a single twin boundary located on one corner; 
this region of material with flipped a, b axes was only ~ 
5% of the total crystal volume. X-ray diffraction was 
used to identify the crystal axes as well as to confirm the 
(essentially) untwinned nature of the crystal. 

The temperature dependence of the penetration depth 
and the superfluid density [p s s = (Ai(0)/Ai(T)) 2 ] for the 
field applied along the a and b axes is shcap. in Fig. [l]. 
A(0) values were taken from Bazov et aZ.BEj (A a (0) = 
1600A and A b (0) =1200 A). It can been seen that n s is es- 
sentially linear down to the lowest temperature measured 
(1.4 K). A fit to the dirty d-wave interpolation formula 
of Hirschfeld and GoldenfeldM [A(T) = aT 2 /(T + T*)} 
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FIG. 1. Temperature dependence of A a and A;, for our un- 
twinned single crystal of YBa2Cua06.95, in zero dc field. The 
superfluid density [p s ,i = (Xi(0)/\i(T)) 2 ] is also shown. 

yields values of T* < IK for both directions. We note 
that the measured anisotropy in dX/dT at low temper- 
atures is only - 3% (d\ a /dT = 4.34 A/K and dX b /dT 
= 4.23 A/K) which is somewhat less than that reported 
by Bonn et aln This may reflect a lower anisotropy in 
the zero temperature values of A in our crystals. We 
note that our values for dp s /dT are within 5% of those 
in Ref. 0. The values of T* for this crystal are among 
the lowest reported to date, showing that our crystal is 
extremely clean. This property is extremely important 
as the Yip-Sauls theory predicts that even small amounts 
of impurities depress the field dependence of A. 

A major problem in measuring the field dependence of 
A, by the method described above, is isolating changes 
in the London penetration depth from other effects 
which also increase the effective penetration depth Xe- 
The most obvious extrinsic contribution is caused by the 
motion of vortices. In the mixed state the additional 
field penetration caused by vortex motion is given by 
Campbell penetration depth A p = B(f>o/ 'p,QH p (k p is the 
pinning force, and <po is the flux quantum). In the limit 
that A p is small, compared to Xl , the contributions add 
in quadrature,E3 X 2 E (H) = X 2 L (H) + X 2 and so 



AAf 



AA, 



2fJ, XLKp 



X v < A; 



(2) 



Using values of k p taken from Ref. |l3j we find that at low 
temperature the vortex contribution would be ~ 25A in 
100 Oe (for H\\ab) \iH c \ was zero and flux was free to en- 
ter the sample. For H less than the actual H c \ the sample 
is in the Meissner state and should, in principle, be free of 
vortices. However, in practice vortices may enter well be- 
fore H c i at sharp corners and surface imperfections and 
so in these regions A^ > A^ • In the measurement config- 
uration H\\c this problem is likely to be intensified as the 
field lines are highly distorted due to the flat plate geom- 
etry of the sample. Another contribution may come from 
weak links at the surface. An applied field may dramat- 
ically suppress the critical current of any weak link thus 
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FIG. 2. Field dependence of A in the a and 6 directions at 
1.4 K. 

causing an increase in A^. Finally in the configuration 
H\\ab part of the supercurrent flows along the c-axis, and 
Xe oc A a (, + t/l\ c (where t and / are the dimensions of 
the sample along the c and a, b axes respectively). As al- 
ready mentioned, because of the relative weak anisotropy 
of YBa2Cu3 6.95 this introduces only a small correction 
to the measured temperature dependence of A for thin 
samples. However, the contribution of A c to the field de- 
pendence of A may not be negligible. From the above it 
is clear that the measured field induced changes in A can 
only set an upper limit for the intrinsic field dependence 
of Al, in all field configurations. 

The field dependence of A for our YBa 2 Cu3 6 .95 crys- 
tal in the a and b directions at our lowest temperature 
(1.4K) is shown in Fig. g. Although AX(H) is approx- 
imately linear with field its magnitude is much smaller 
than the Yip-Sauls prediction. The slope of the straight 
line fitting through that data is 20 ± 5 x 10~ 3 A/Oe (i.e., 
~ 2A in 100 Oe), and within the error is the same in both 
directions. This is more than a factor 2 times lower than 
the Yip-Sauls prediction. No hysteresis was seen within 
the scatter of the data. 

Data for the changes in the average in-plane penetra- 
tion depth (X a b) m the third field configuration (H\\c) 
are shown in the upper panel of Fig. [| AX a b(H) is lin- 
ear up to the maximum field measured [a decrease in 
slope was seen above ~ 20 Oe (not shown)]. Again no 
hysteresis was observed at these low field values. For 
this flat sample, in this geometry, the field at the surface 
of the sample is clearly very different from the applied 
field. We have estimated the field at the edge of the 
sample by multiplying the applied field by a demagnetiz- 
ing factor [1/(1-N)], derived from an inscribed ellipsoid 
approximation. This 'demagnetizing' factor was within 
20% of that estimated by comparing the known volume 
of the sample to the measured apparent shielding vol- 
ume (as measured in a zero field cooled dc magnetization 
measurement). The actual field on the edge of the sam- 
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FIG. 3. Top: Field dependence of X a b with the field applied 
along c at 1.4K. Bottom left: The evolution of dX a b/dH with 
temperature. Theory curves are the average slope AX(H)/H 
with H = 100 Oe and 300 Oe from Ref. |. Bottom Right: 
Fraunhofer like oscillations in A(_ff) in another sample. The 
data are highly reproducible (5 field sweeps are shown). 

pie will not be constant over the whole width, however 
this 'demagnetizing' factor is a rough estimate of the av- 
erage field. Using 1/(1-N)=36, we estimate the slope 
dX a b/dH = 6.3 x 10~ 3 A/Oe, which is ~ 3 times smaller 
than for the H\\a, b measurements and more than 7 times 
smaller than the Yip-Sauls prediction. We note that with 
this demagnetizing factor that the change in slope at H~ 
20 Oe corresponds roughly to H c \ in this direction. 

The lower left panel of Fig. || shows the evolution of 
the slope, dX a b/dH as temperature is increased. There 
is almost no change in dX a b/dH up to 50K, in sharp con- 
trast to the predictions of the Yip-Sauls theory. As T c is 
approached dX a b/dH increases substantially. There was 
no indication of the predicted change from AX(H) oc H 
to AX(H) ex H 2 as temperature was increased. We con- 
clude from these observations that the small field depen- 
dence that we see does not originate from the Yip-Sauls 
effect. If we interpret AX(H) as originating from vor- 
tex motion, the above results would imply that at low 
temperature k p is almost constant in our samples. 

A close examination of X a b(H) in Fig. || shows that the 
apparent noise on top of the linear dependence is repro- 
ducible with a pattern which is approximately symmetric 
about H = 0. In fact these modulations in X a b are quite 
reproducible and become more pronounced as tempera- 
ture is increased. In another sample we observed much 
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FIG. 4. Field dependence of X a for a second sample, at 
several fixed temperatures (indicated on figure). H was \\b. 
The inset shows the temperature dependence of p a>a and p a ,b 
for the same sample. 



bigger modulations which resembled closely the type of 
Fraunhofer pattern typical of the field dependence of the 
critical current in weak links (see the lower right panel in 
Fig. H). We conclude that these effects indicate the pres- 
ence of weak links on the edge of the crystal, the critical 
current of which is modulated by the static dc field. We 
note that this structure is only observable because of the 
extremely low modulation fields used in our experiment. 

Finally in Fig. |] we show data for a second sample 
of optimally doped, untwinncd YBa2Cu3 06.95 (T c =93.6 
K). p s (T) (inset Fig. ^) for this sample is linear down to 
only ^4K, showing that it is dirtier than the first sample 
(cf., Fig. [j]). AX(H) however is larger by around a factor 
2 and is linear up to a field of approximately 50 Oe. 
Again no hysteresis was observed over the field range 
in the figure (H <100 Oe). This linear slope increases 
slightly with increasing temperature in sharp contrast to 
the predictions of the Yip-Sauls theory. From this fact we 
conclude that this linear contribution is again of extrinsic 
origin. The larger magnitude of AA(iJ) probably reflects 
more flux leaking into this sample. 

The weak temperature dependence of dX a b / dH shown 
here is consistent with the results of Maeda et alB for 
their samples of YBa2Cu30g.95 and Bi2Sr2CaCu208- In 
their work the lowest temperature at which experiments 
were performed was 10K, and little change was found 
in dX/dH as temperature was increased. However these 
authors report much bigger changes in X(H) than are 
found here or are indeed predicted by the Yip-Sauls the- 
ory. The large size of these effects and the fact that they 
are not damped out by increasing temperature leads us 
to conclude that they do not originate from the Yip-Sauls 
effect, and probably arise from vortex motion. 

In summary, we have presented measurements of the 



field dependence of A with a resolution around a fac- 
tor 100 better than previous experiments, on very high 
purity single crystals of YBa2Cu3 6.95, at temperatures 
down to T/T c =0.015. In our cleanest sample we observe 
a very weak field dependence in all three crystallographic 
directions. The overall magnitude of this change in A is 
at least two times smaller ( or 7 times smaller in the H\\c 
configuration) than our estimates of what we would ex- 
pect to see if the Yip-Sauls theory was valid. A second 
sample shows a larger effect even though its purity level 
is lower. In all cases however the observed field depen- 
dence of A is not diminished by increasing temperature 
and so we conclude that it does not arise from the Yip- 
Sauls effect. The sample dependent AX(H) we observe 
most likely originates from vortex motion. 

The large body of evidence collected thus far shows 
fairly convincingly that YBa2Cu3 06.95 is a d-wave su- 
perconductor and so other possible reasons must be con- 
sidered for the apparent absence of the Yip-Sauls effect 
below H c \. One possibility could be the effect of non- 
local electrodynamics. It is usually assumed that as the 
in-plane penetration depth (Ao ^1600A) is much greater 
than the in-plane coherence length (£o ~ 16A)_then local 
electrodynamics are sufficient. Kosztin et al&3 however 
have pointed out that this argument does not hold for a 
rf-wave superconductor, as £ is angle dependent. At the 
nodes the energy gap A(0) cx £(6*) -1 goes to zero and 
£(0) diverges violating the local condition. As the region 
near to the nodes is exactly where the low temperature 
excitations are formed, this observation can be expected 
to have serious consequences for hoth the temperature 
and field dependence of A. Li et al.lB have recently made 
calculations of the AA(-ff) which include both the non- 
local and non-linear effects. They conclude that the non- 
local effects indeed drastically reduce AX(H) at fields 
< H c \, rendering the effect essentially unobservable in 
the Meissner state. Although these calculations strictly 
only apply in the H\\c configuration (when the direction 
of field penetration is in the same plane as the super- 
currents), the relatively low anisotropy of YBa2Cu30g.95 
means that a similaweduction may also occur in the H\\ 
plane configuration.!^ 
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